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We discuss and model an experiment to study quasi-condensate growth on an atom chip. In 
particular, we consider the addition of a deep dimple to the weak harmonic trap confining an 
ultracold one-dimensional atomic Bose gas, below or close to the characteristic temperature for quasi- 
condensate formation. The subsequent dynamics depends critically on both the initial conditions, 
and the form of the perturbing potential. In general, the dynamics features a combination of 
shock-wave propagation in the quasi-condensate, and quasi-condensate growth from the surrounding 
thermal cloud. 
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I. INTRODUCTION 

The recent achievement of effectively one-dimensional 
ultracold atomic samples on microfabricated surfaces, 
known as 'atom chips' lj opens up the way for a num- 
ber of applications, such as precision intcrferometric 
measurements H H H and quantum computation 
H H 01 One-dimensional (ID) geometries, in which 
the transverse confinement exceeds all other relevant en- 
ergy scales in the system, offer improved atomic guid- 
ance and device miniaturization d EJ- This comes 
at the expense of increased phase fluctuations, which 
tend to destroy the coherence of the ensemble. Var- 
ious equilibrium studies have already been performed 
both experimentally O JigjP J3and theoretically 
El 13 Hi 13 ES l2fl ET 121 EST 13 to understand 
in particular this latter effect. However, a deeper under- 
standing of the intrinsic dynamics of such systems is still 
lacking, including the issue of quasi-condensate forma- 
tion, and the role of phase fluctuations on its growth dy- 
namics. Although condensate formation has been stud- 
ied in 3D geometries |3 EI El El El H3 , such studies 
are only now possible in ID. In this paper, we discuss 
a realistic experiment which captures the dynamics of 
quasi-condensate formation and relaxation into a per- 
turbed potential, by adding a dimple microtrap on an 
atom chip containing a gas of ID ultracold Bose atoms. 

The issue of quasi-condensate formation may be re- 
lated to a recent 3D condensate growth experiment E3 , 
which to date remains largely not understood. In par- 
ticular, this experiment revealed, under conditions of 
slow cooling, an unexpected slow linear initial conden- 
sate growth. As these features were observed close to the 
transition point, the authors suggested this might be due 
to the enhanced phase fluctuations in this region, lead- 
ing to the formation of a quasi-condensate that preceded 
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the usual condensate growth. The dimensionality of our 
envisaged ID experiment ensures the system remains in 
the regime of large phase fluctuations throughout its en- 
tire evolution. Moreover, under conditions of slow cooling 
from above the transition point, our results point towards 
a slow initial linear-like growth which is consistent with 
bosonic enhancement. However, a detailed comparison to 
the above experiment would require a full 3D calculation, 
and therefore lies beyond the scope of this paper. 

The experimental scenario considered here is the fol- 
lowing: An atom chip is first loaded with a certain 
number of atoms, at densities low enough to be in the 
ID regime, and cooled to a prescribed low temperature. 
The approximately harmonic longitudinal confinement of 
the atoms is then perturbed by the addition of a deep 
dimple microtrap, of variable width, as shown in Fig. 
1. Such a technique of local phase-space compression 
[3ll |3| has been used to force a system of 23 Na atoms 
into the quantum-degenerate regime in a reversible man- 
ner in 3D [3|- A related approach has been used to 
create quantum-degenerate 133 Cs in 3D [313]' an d in 
2D [35|]. In ID, such a dimple can be created by optical 
and magnetic traps, although the optimum technique on 
an atom chip appears to be the application of electric 
fields to a separate electrode, as discussed in Appendix 
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FIG. I: Initial (brown lines) and final (black) longitudinal 
harmonic potentials (solid) and total densities (dashed) for a 
system of 2500 87 Rb atoms at a temperature T = 50nK for 
(a) a wide dimple w = 41 z , and (b) a tight dimple w = 0.5/ z , 
with zd = 0, Vo = —2000huj z . Here, l z — ^h/mu z is the 
longitudinal harmonic oscillator length. 
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A. We study the resulting growth dynamics focusing on 
the regime of a very deep dimple, whose potential greatly 
exceeds all other energies in the system , as this leads 
to the most interesting non-equilibrium regime. If there 
is already a quasi-condensate present in the original trap, 
the addition of the dimple is found to lead to shock-wave 
formation and subsequent relaxation of the perturbed 
quasi-condensate in the dimple. In the opposite case of 
an initially incoherent ultracold atomic sample, we ob- 
serve direct quasi-condensate growth from the surround- 
ing thermal cloud. The nature of the dynamics addi- 
tionally depends on the width of the dimple compared to 
the size of the atomic cloud, and on the dimple location. 
The distinct regimes of quasi-condensate dynamics can 
be probed by controlling the initial temperature, atom 
number and perturbing potential. 

This paper is structured as follows: Sec. II discusses 
the methodology used to model the envisaged experi- 
ment. The dynamics following the addition of the dimple 
in the limit of low temperature is then discussed in Sec. 
Ill, paying particular attention to shock-wave formation 
(Sec. Ill A), its effect on the dynamics of the density at 
the trap centre (Sec. Ill B), and the effect of the width 
of the dimple compared to the size of the atomic cloud 
(Sec. Ill C). The effect of temperature is discussed in 
Sec. IV for both wide and tight dimples, with the latter 
case highlighting the interplay between spontaneous and 
stimulated growth in the dimple. Sec. V summarizes the 
examined regimes by means of a suitable graph, and Sec. 
VI features a brief summary. The discussion of experi- 
mental techniques which can be used to create the dimple 
trap on an atom chip has been deferred to Appendix A. 



II. METHODOLOGY 

To model the proposed experiment, we first prepare 
the desired initial state on the atom chip stochastically. 
This is achieved by means of the Langevin equation 
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describing the dynamics of the order parameter 
$(z,t) [HSUII. Here KxtO) = muj 2 z 2 /2 denotes the 
longitudinal harmonic confinement, and g = 2haujj_ is 
the one-dimensional coupling constant obtained by aver- 
aging over transverse gaussian wavefunctions, where a is 
the three-dimensional scattering length and uj± the trap 
frequency in the transverse directions [40l | . The chemical 
potential, /it, determines, for a given initial trap potential 
and temperature, the total atom number in the system. 
The atom chip trap is pumped from a thermal reservoir 
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In accordance with the fluctuation-dissipation theorem, 
there is an associated gaussian noise contribution r](z, t) 
obeying (rj*(z,t)rj{z',t')) = (ih 2 /2)T, K (z)S(z~z')S(t-t'), 
where the brackets denote averaging over the realizations 
of the noise. The dependence of these quantities on the 
Keldysh self-energy TiYj K (z) ensures that the trapped gas 
relaxes to the correct thermal equilibrium, as addition- 
ally verified by direct comparison to the modified Popov 
theory discussed elsewhere 0, El 113 ■ F° r simplic- 
ity, we choose in this work both the longitudinal and 
the transverse atom chip confinement to be harmonic, 
with respective trap frequencies uj z = 2tt x 5 Hz and 
uj± = 2tt x 5000 Hz, and present explicitly results for 
87 Rb, with a — 5.32nm. Unless otherwise specified, re- 
sults will be given in dimensionless harmonic oscillator 
units of the original trap, such that lengths are s caled to 
the longitudinal harmonic oscillator length l z — ^Jh/mui z 
and energies to huj z . 

Once the system has fully relaxed to the desired equi- 
librium, we suddenly turn on a gaussian dimple potential 
and simultaneously remove the coupling to the heat bath. 
Switching off driving and dissipation terms, i.e., taking 
R(z, t) = ■q(z 1 1) = 0, ensures that the total atom number 
on the atom chip remains fixed, as in the envisaged exper- 
iment, and thus only intrinsic dynamics of the system are 
taken into account. The dimple potential added is given 
by Vd = Vbe~( z ~ ZD ' l 2w where Vb is the dimple depth, 
w its width, and zd its location. A sudden introduction 
of the dimple potential is chosen, since this generates a 
highly non-equilibrium situation, which gives rise to the 
most interesting dynamics. For the chosen parameters, 
this corresponds to a trap turn-on timescale of 6/xs, which 
is an experimentally realistic timescale. The subsequent 
analysis is based on averaging over many different initial 
realizations, which ensures that both density and phase 
fluctuations are accurately included. This is equivalent to 
averaging over a large number of independent experimen- 
tal realizations with a variable initial phase, as typi cally 
done when investigating growth dynamics [25t [29. |27| . 
In our discussion, we assume the dimple is turned on at 
t = 0. 



III. LOW TEMPERATURE LIMIT 

The sudden addition of a deep gaussian dimple on a 
harmonically confined pure condensate leads to a large 
atomic flux towards the center of the dimple, and thus 
to the development of a large local density gradient at 
symmetric points about the dimple centre, which, in 
turn, leads to the formation of two counter-propagating 
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shock wavefronts. The dynamics of the initial shock 
wave formation stage for a pure coherent condensate has 
been discussed under various related conditions in Refs. 
|4ll.l43 .l43j. In addition, shock waves have already been 
observed experimentally in the context of rapid poten- 
tial perturbations related to the formation of dark soli- 
tons 01 an d vortex lattices [4^. In this paper, we dis- 
cuss both short and long term dynamics in the dimple, 
and identify the distinct growth dynamics induced by 
the addition of the dimple. In particular, we highlight 
the competing effects of temperature, chemical potential 
and dimple width on these processes. 

We start our analysis by considering the most general 
scenario of experimental relevance. The dimple of Fig. 
1(a) is added onto a quasi-condensate with about 3,400 
atoms. This particular dimple has been chosen, as it can 
be easily generated with existing atom chips. A detailed 
discussion of how this can be achieved can be found in 
Appendix A. 

In brief, we find the subsequent dynamics to be es- 
sentially controlled by two parameters: (i) Firstly, the 
dynamics depends critically on the amount of quasi- 
condensation present in the initial system. This is de- 
termined by the ratio of T/T c , where T c is the effective 
ID 'transition' temperature, (ii) Secondly, the dynam- 
ics is sensitive to the relation of the spatial extent of the 
atomic cloud in the original trap compared to the effec- 
tive width of the perturbing dimple potential. 



A. Non-equilibrium Dynamics in the Dimple 

For simplicity, we focus initially on the regime T <C T c , 
for which most of the atoms are in the quasi-condensate. 
Typical density snapshots of the growth dynamics in this 
limit are shown by the black lines in Figs. 2 (a)-(i). In 
our analysis, the 'transition' temperature T c in the pres- 
ence of interactions is determined numerically from our 
modified Popov theory |l6|, Il7l llSL Il9| . 

Once the dimple is turned on, the central density in- 
creases rapidly, until the instability towards shock-wave 
formation is reached |4lJ, |42j, |43j. The ensuing shock 
waves propagate towards the dimple edge and subse- 
quently reflect back into the dimple centre. Individual 
profiles can be simulated by solving the Gross-Pitaevskii 
equation for a fully coherent condensate, yielding den- 
sity profiles with large density variations, as shown by 
the green curves in the insets of Figs. 2(d)-(h). How- 
ever, a detailed understanding of growth dynamics which 
enables a straightforward comparison to experiments re- 
quires the addition of random initial phases and subse- 
quent averaging over such initial configurations. Such 
averaged profiles, obtained within our analysis, reveal 
the initial formation of a large central peak, followed 
by its break-up into two smaller peaks, moving in oppo- 
site directions towards the dimple edges, with the central 
peak re-emerging after a characteristic time determined 
by the dimple potential. The latter peak further splits 
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FIG. 2: (color online) Averaged snapshots of total density in 
the initial harmonic (a) and perturbed ((b)-(i)) trap, taken 
at equal time intervals, At = 0.073uf^~ , for a gas of 3400 
87 Rb atoms. Profiles are shown at T = 50nK < T c (black 
lines) and T = 200nK ~ T c (brown), (i) Indistinguishable 
equilibrium profiles due to large dimple depth. Insets com- 
pare single-run results based on the Gross-Pitaevskii equation 
(green) to averaged stochastic profiles at T = 50nK. (j) Cor- 
responding oscillations in averaged density at the trap centre. 
Inset highlights the first two oscillations, showing explicitly 
the times of the plotted snapshots (a)-(h), and comparing to 
results of Gross-Pitaevskii (green). Dimple parameters as in 
Fig. 1(a). 



into counter-propagating lower peaks, and so forth, un- 
til the gas equilibrates to the profile shown in Fig. 2(i). 
Such shock-wave dynamics would be largely suppressed 
if the dimple was much shallower, or if it were introduced 
on a much slower timescale. 

The long-term system dynamics is portrayed in the 
evolution of the central density of the gas shown by the 
black lines in Fig. 2(j). It features oscillations, on top 
of a growth curve. This complicated dynamics is a re- 
sult of a number of competing processes occuring simul- 
taneously. Firstly, the observed oscillations are a direct 
consequence of averaged shock-wave propagation in the 
quasi-condensate. In fact, the initial oscillatory dynamics 
in the low temperature limit resembles a suitably aver- 
aged single-run of the Gross-Pitaevskii equation, shown 
by the green lines in the inset to Fig. 2(j). The 'growth' 
part of the dynamics arises as a result of spatial com- 
pression of the trap, and features various contributions: 
In addition to quasi-condensate compression, there is an 



4 



additional process whereby all thermal atoms located in 
the dimple area fall into the perturbed trap by joining 
the quasi-condensate. This is ensured by the depth of 
the applied dimple, which largely exceeds all other rel- 
evant energies of the system. Moreover, thermal atoms 
located in the tails of the initial atomic cloud (as well as 
any quasi-condensate atoms located outside of the region 
of the dimple) are pulled into the trap centre, continu- 
ously interacting with the propagating quasi-condensate 
shock wave, and eventually leading to additional quasi- 
condensate growth. In the remaining part of the paper we 
discuss these competing processes in more detail by iden- 
tifying suitable experimentally realistic limiting regimes. 



B. Oscillation Frequencies 

Firstly, we comment on the observed oscillations. The 
addition of the dimple corresponds essentially to an in- 
stantaneous local compression of the central region of the 
trap. For an easy visualization which captures the main 
dynamics discusssed here, let us consider the simplified 
case of a harmonic dimple of width comparable to the 
system size, such that essentially all of the dynamics is 
contained within a suitably defined 'dimple region'. In 
this case, the effect induced by the addition of the dim- 
ple is similar to that of an instantaneous increase in the 
harmonic frequency. This leads to excitation of the low- 
est compressional mode, which, for a pure condensate in 
the ID mean field regime, has been predicted to occur 
at a frequency of V3wd Ha,|43, as already observed ex- 
perimentally |48j . and additionally verified numerically 
in our simulations. 

In the experimentally more realistic gaussian dimples 
considered here, we can define an effective harmonic 
frequency uod — VVo/w at the central dimple region. 
Although the initial oscillation in the gaussian dimple 
occurs very close to the predicted frequency, the fre- 
quency shifts towards lower values, as soon as the quasi- 
condensate no longer feels a purely harmonic confine- 
ment, i.e., for system sizes larger than the effective dimple 
width. This shift becomes more pronounced in time as 
the quasi-condensate grows in the dimple. In addition, 
the observed oscillations are damped due to various ef- 
fects. Firstly, unlike for a purely coherent system, the 
averaging performed over random initial phases in the 
presence of fluctuations leads to dissipation, precisely as 
would be the case in averaged condensate growth ex- 
periments in the presence of a non-negligible thermal 
cloud. Furthermore, the relative motion between the 
shock waves and the (mostly thermal) atoms entering 
the dimple region from the edges is expected to create 
an additional channel for dissipation. Finally, the addi- 
tion of a gaussian dimple leads to excitations of multiple 
frequencies, which lead to beating effects and an acceler- 
ated decay of the oscillation amplitude. 
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FIG. 3: (color online) (a) Initial density for a gas of N x 960 
atoms at T = 25nK and fi = 100hu> z , such that (i/ksT sa 
1 (solid blue lines), versus corresponding zero-temperature 
Thomas- Fermi profiles (dashed green lines). The harmonic 
potential is also shown for easier visualization (solid black), 
(b) Initial (solid black) versus final (solid brown) potentials for 
the two dimples of Fig. 1. Corresponding dimple potentials 
Vd are shown by the dashed red lines. 



C. Effect of Dimple Width 

The above example featured a combination of shock- 
wave propagation and growth. These effects can be iso- 
lated by changing the form of the perturbation. Shock 
waves arise in the quasi-condensate already present prior 
to the addition of the dimple. On the other hand, the 
observed increase in the central density is the result of 
both local compression of the atoms in the dimple re- 
gion, and of direct growth arising from (mainly thermal) 
atoms moving into the central region from the trap edge. 
As a result, oscillations dominate when the dimple width 
is comparable to the quasi-condensate system size, and 
for relatively small atom numbers, such that growth is 
minimized. On the other hand, the growth features are 
dominant in very tight traps for which the dimple width 
is much smaller than the system size. 

These two contrasting regimes for the two opposite lim- 
its of wide and tight dimples, with respect to the effective 
system size, are shown in Fig. 3 for a reduced number 
of N w 960 atoms for computational convenience. To 
facilitate an easy comparison between these two limits 
we keep the dimple depth fixed. To give a simple vi- 
sualization of this distinction, we note that, for the low 
temperature case considered here, the zero-temperature 
Thomas- Fermi profile (dashed green line in Fig. 3(a)) is 
a reasonable first approximation to the system density 
prior to the addition of the dimple (solid blue line). In 
the chosen harmonic oscillator units, the zero temper- 
ature Thomas- Fermi radius is defined by Rtf = V^fJ-, 
where fj, is the chemical potential of the system. 

An effective dimple width, Hp can also be approxi- 
mately defined as the point at which the dimple depth 
falls to 0.01 of its maximum value Vo, implying Rd w 3w. 
Fig. 3 plots both initial densities and potentials (top) 
and perturbed potentials (bottom) for the two opposite 
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FIG. 4: (color online) (a)-(b) Dynamics of dimple central 
density for N ~ 960 atoms upon adding a broad dimple with 
w — 4l z (black) or a narrow one with w = 0.5l z (brown), as 
in Fig. 3. Time plotted in terms of (a) original trap timescale 
uiz 1 and (b) dimple timescale cj^ 1 , with lod = 11.2oj z (black), 
or 89.4tj z (brown). Here T = 25nK and fi = 100hu z , with 
fi/ksT ~ 1. Insets: (a) Same plot with densities approxi- 
mately scaled to final equilibrium values iieq{0), and (b) Ini- 
tial regime highlighting the overlapping density oscillations. 



regimes of Rd ~ Rtf (left) and Rd <C Rtf (right). 
Note that the equilibrium densities in both dimple traps 
are essentially given by the Thomas-Fermi profiles in 
these traps. This is true for all atom numbers consid- 
ered in this paper (N < 3500), for which the chosen large 
dimple depth ensures that all atoms can be accomodated 
within the central approximately harmonic dimple region 
(see also the density profiles in Fig. 1). 

The system dynamics in these two opposite regimes 
is portrayed by the evolution of the central density in 
Fig. 4. This evolution is plotted both in terms of the 
original trap timescale uj~ 1 in Fig. 4(a), and in terms 
of the effective dimple timescale uij-, 1 in Fig. 4(b), as 
these reveal different dynamical features that we wish to 
comment on. 

The actual relaxation timescale, and long-term dynam- 
ics in these opposing regimes are best compared when 
both curves are plotted in terms of the characteristic 
timescale cjJ 1 of the initial harmonic trap. Fig. 4(a) 
thus shows that the central density in the tight dimple 
increases much more significantly, due to the enhanced 
spatial compression. A comparison of the growth rate 
in the different dimple traps is facilitated in the inset to 
Fig. 4(a), by plotting the dynamics of the peak densi- 
ties, approximately scaled to their respective equilibrium 



values. Although equilibration timescales appear to be 
comparable, the initial dynamics is slower in the tight 
trap, presumably due to the reduced overlap between the 
states in the initial and final traps. 

The initial dynamics, including the shock-wave- 
induced oscillations are more appropriately investigated, 
when the same results are plotted in terms of the effective 
dimple timescale u>d , as in Fig. 4(b). This is because 
the oscillations are actually fixed by the frequency u>d 
of the perturbed trap (multiplied by the factor of V3), 
which varies with dimple width w. When plotted in this 
manner, the initial dynamics feature a striking similar- 
ity in both amplitude and frequency of the respective 
central density oscillations, as shown in the inset. This 
is due to the fact that the initial dynamics is set by the 
quasi-condensate which is, in both cases, present over the 
entire extent of the dimple region. However, the subse- 
quent growth is noticeably different, due to the enhanced 
compression in the tight dimple, leading to a much en- 
hanced growth in the central density. 

IV. EFFECT OF TEMPERATURE 

All earlier discussion focused on the regime T <C T c . 
We now discuss the effect of changing the initial temper- 
ature of the system, with respect to T c . 



A. Wide Dimple 

The density snapshots shown in Fig. 2 correspond to 
fairly 'typical' experimental regimes, for which Rtf{0) 
is approximately equal to a few Rd, such that the ad- 
dition of the dimple is accompanied by a combination 
of oscillations and growth. However, as the temperature 
of the initial sample increases towards T c at fixed atom 
number, the propagating secondary density peaks in the 
averaged profiles become washed out, as shown by the 
brown lines in Fig. 2(a)-(h). This is due to the enhanced 
fluctuations in the initial state prior to the addition of 
the dimple. Nonetheless, the final averaged profile is al- 
most independent of temperature, due to the large dim- 
ple depth, which overshadows all other energies of the 
system. 

By comparing these two cases, we note that the am- 
plitude of the initial oscillation is controlled by the ra- 
tio of the chemical potential [i to the thermal energy 
ksT. Under conditions of strong condensation, such that 
fi/ksT > 1, the initial oscillation amplitude is large and 
can even exceed the equilibrium value. In the opposite 
regime, its amplitude is largely suppressed. In the limit 
of a sufficiently small atom number and a relatively broad 
dimple, the addition of the dimple perturbs the central 
density only mildly. Thus, the additional constraint of 
extremely low temperature (fcsT <C fJ.), leads essentially 
only to the appearance of oscillations without substantial 
growth. 
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B. Tight Dimple 

The tight dimple, whose width is much smaller than 
the effective system size, creates a natural temporal sepa- 
ration between the initial shock- wave-induced oscillations 
and the subsequent growth dynamics. These competing 
effects exhibit different characteristics, depending on the 
ratio T/T c , as shown in Fig. 5(a). This ratio can be 
controlled either by changing the temperature at con- 
stant total atom number, as discussed for Fig. 2 above, 
or by varying the atom number at fixed temperature. 
In this section, we choose to discuss the latter, as this 
enables a further distinction between spontaneous and 
stimulated bosonic growth. Note that, for the small atom 
numbers 400 < N < 1000 considered here, the condition 
Rtf(0) Rd is always satisfied for the dimple with 
w = 0.5l z . 

Regarding the initial oscillations in the central density, 
we note that, as the atom number decreases, so does also 
the ratio of /x/fc^T determining the amplitude of the first 
oscillation. This leads to very strong suppresion when 
fi/ksT <C 1, as shown by the bottom curve in the inset to 
Fig. 5(a). After quenching of these initial oscillations, we 
observe a 'secondary growth' dynamical phase. For small 
atom numbers, corresponding to large T/T c , we observe 
a very slow initial growth in the central density. This 
will be shown to be consistent with spontaneous quasi- 
condensate growth from the thermal cloud. Increasing 
the atom number leads to a lower value of T/T c , and to 
a significant quasi-condensate fraction in the initial trap, 
which yields a faster initial growth rate in the dimple. 

To understand this in more detail, Fig. 5(b)-(c) plots 
the integrated quasi-condensate atom number in the dim- 
ple as a function of time for the cases (b) T w T c , and 
(c) T -C T c . The quasi-condensate atom number is deter- 
mined from a combination of local quadratic and quar- 
tic correlations of the order parameter &(z,t), as will 
be discussed in more detail elsewhere. The evolution is 
in both cases fitted by two distinct growth models, as 
in [25| . The first model includes both spontaneous and 
stimulated bosonic growth, and is fitted by dashed brown 
lines. It obeys 
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N Q (t) = N» 

where Nq the total quasi-condensate atom number in the 
dimple, the superscripts (i) and (f) denote respectively 
initial and final values, and 7 corresponds to the initial 
growth rate. The value 6 = 2/3 assumes a growth rate 

linear in the difference of chemical potential of the ther- 

2/3 

mal cloud and the quasi-condensate, for which /1 oc N . 
The second model, fitted by the solid green lines, de- 
scribes simple exponential relaxation to equilibrium via 



N (t) = iV (/) (1 - e-f*) 



(4) 




FIG. 5: (color online) (a) Central density dynamics in tight 
dimple (it) = 0.5i z ) with increasing atom number (from bot- 
tom to top) N w 420, 490, 620 and 960, with R D < R TF in 
all cases. The temperature is fixed to T — 25nK, and corre- 
sponding chemical potentials are fi/%w z = 5,30,60,100. In- 
set: Comparison of initial oscillatory dynamics, (b)-(c) Quasi- 
condensate growth in the dimple for (b) N « 420 (T ~ T c ), 
and (c) N « 960 (T < T c ). Each computed curve (black) 
is fitted with the growth curve of Eq. (3) (dashed brown), 
and the relaxation curve of Eq. (4) (solid green). The initial 
number of condensate atoms in the dimple is bigger in (c), 
due to the larger chemical potential. 



The limit T w T c shown in Fig. 5(b) is well described 
by the growth model of Eq. (3), indicating spontaneous 



and subsequent stimulated quasi-condensate growth in 
the dimple, analogous to the first condensate growth 
studies in 3D systems. Although the exponential re- 
laxation model fails to describe the entire dynamics in 
this case, it nonetheless works well if fitted after a cer- 
tain 'initiation' time [2^|. In the opposite regime of 
r«r c shown in Fig. 5(c), the initial dimple dynamics 
is mainly governed by re-equilibration of the perturbed 
quasi-condensate in the combined trap. This obeys the 
exponential relaxation model well, indicating that the 
dynamics is dominated by stimulated growth. 



The intermediate temperature regime features inter- 
esting, but complicated dynamics, as the effects of spon- 
taneous and stimulated growth compete with each other. 
The above distinction between spontaneous and stimu- 
lated growth becomes less pronounced with increasing 
atom number, with the initial growth dynamics being 
well-modeled by an exponential relaxation curve even for 
the tight dimple considered here, when N is a few 1000 
atoms. 
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V. SUMMARY OF DISTINCT REGIMES 

The effect of adding a dimple trap to the centre of a 
weaker harmonic trap containing a ID quasi-condensate 
was considered in the limit when the dimple depth greatly 
exceeds all other relevant energies in the system, i.e., 
Vb 3> Mi ksT, huj z . The presented analysis was restricted 
to the typical experimental regime satisfying the follow- 
ing length scale separation 

1±<Z< (A d s, Id) « lz « (Rd, Rtf(0)) (5) 

where l± = ^h/mujj_ is the transverse harmonic oscil- 
lator length, £ = fi/^JA^m is the healing length of the 

system, XdB = \J 2irh 2 /rnksT is the thermal dc Broglie 

wavelength, Id — \/h/mu)D is the effective harmonic 
dimple oscillator length, Rtf(0) is the zero-temperature 
Thomas-Fermi radius in the original harmonic trap, and 
Rd is the effective dimple width. Within this regime, 
the dynamics in the dimple displays an interesting inter- 
play between shock-wave propagation in the perturbed 
quasi-condensate, and direct quasi-condensate growth. 

The following important conclusions were reached 
about the accessible dynamical regimes in such systems: 

• The initial dynamics for T < T c is dominated by 
shock-wave propagation, leading to large oscilla- 
tions in the central density. The initial amplitude 
of such oscillations is controlled by the ratio of 
(n/ksT), yielding large oscillations when this is of 
order unity or larger. 

• The long-term system dynamics is dominated by 
quasi-condensate compression and growth in the 
dimple, with this effect largely pronounced when 
the system size greatly exceeds the dimple width. 

The different accessible regimes, along with the asso- 
ciated relevant parameters for (/x/fcsT) and Rtf(0)/Rd 
are summarized by the evolution of the central density 
shown in Fig. 6. In this figure, each horizontal line (a)-(c) 
corresponds to a particular configuration of fixed initial 
atom number and dimple width. Left images (i) show 
the case r<T c , while right images (ii) display the high- 
temperature limit T « T c . The images shown portray, 
from top to bottom, the following cases: 

(a) A relatively large atom number N = 3400, such 
that Rtf{0) is slightly larger than Rd- This 
leads to noticeable growth, and additionally fea- 
tures pronounced and 'long-lived' oscillations when 
fi/k B T w 1. 

(b) A sufficiently small atom number N = 540 and 
comparatively wide dimple, Rtf{0) ~ Rd, for 
which growth is suppressed and oscillations are the 
dominant feature in the low temperature limit. 




t(0 



FIG. 6: Visualization of distinct regimes accessible in ul- 
tracold ID Bose gases, upon addition of a deep dimple 
(Vo = — 2000?kJ z ) at the centre of a weaker harmonic trap 
(i.e. zd =0), as manifested by the dynamics in the central 
density oscillations. From top to bottom, iV fid (a) 3400, (b)- 
(c) 540, and w = (a)-(b) 4l z , and (c) 0.5/z. Left images (i) 
correspond to T <C T c , with right images (ii) showing the 
opposite regime T fis T c . Each figure displays the correspond- 
ing approximate values of a = h/UbT, and /3 = Rtf(0)/Rd, 
determining the system dynamics. 



(c) A very tight dimple Rd <C Rtf{0), featuring 'long- 
term' pronounced growth, after damping of any ini- 
tial oscillations occuring. Note that the sponta- 
neous bosonic growth contribution is pronounced 
here due to the small number of atoms. 

Having identified the importance of the relation of the 
size of the atomic cloud compared to the dimple width, 
and of the chemical potential compared to the tempera- 
ture of the system, we can now 'tailor' our experimental 
conditions to produce the desired dynamics. 

Note that all dimples were considered to have the same 
depth, which exceeds all other relevant energy scales in 
the system. Thus, for the atom numbers considered here, 
the final densities are all described by Thomas- Fermi pro- 
files in the dimple, and the dimple width is the relevant 
parameter determining the ensuing ID dynamics. 

The case of an off-centered dimple, with z ^ 0, yields 
qualitatively similar behaviour, that will be discussed in 
more detail elsewhere. When the dimple falls within 
the initial quasi-condensate size, shock wave formation 
plays a key role in the subsequent dynamics. However, 
in the opposite limit of a dimple located outside of the 
quasi-condensate, the oscillations are largely suppressed, 
and quasi-condensate relaxation competes with quasi- 
condensate growth. The whole process is further modi- 
fied by the asymmetry in the initial density distribution 
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with respect to the dimple centre, an effect which be- 
comes particularly pronounced for shallow dimples. 



VI. CONCLUSIONS 



Starting from a ID potential we can create a dimple 
on an atom chip by using (i) magnetic fields, (ii) electric 
fields and (iii) dipole potentials. In the following dis- 
cussion of creating a dimple we neglect the longitudinal 
confinement of the original trap. 



In conclusion, we studied one-dimensional quasi- 
condensate growth dynamics in a dimple microtrap cre- 
ated on top of the harmonic confinement of an atom chip. 
Rich novel dynamics were observed in the case of a deep 
dimple, whose depth exceeds all other energies of the sys- 
tem, displaying an interesting interplay between shock- 
wave propagation in the perturbed quasi-condensate, and 
direct quasi-condensate growth. The optimum observa- 
tion of shock-wave dynamics, without the added signifi- 
cant growth of the central density, requires low temper- 
atures T <^iT Cl a comparatively large chemical potential 
fi > ksT, and a dimple of width comparable to the sys- 
tem size. Understanding these processes in more detail is 
expected to contribute to fundamental issues in the dy- 
namics of degenerate one-dimensional Bose gases, with 
potential applications in atom lasers and atom interfer- 
ometers. 
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APPENDIX A: EXPERIMENTAL CREATION OF 
ID DIMPLE POTENTIALS 

This Appendix discusses various experimental tech- 
niques which can be used to create a dimple potential, in 
a quasi-lD atomic gas on an atom chip [l|, Eii lMI > m or- 
der to observe the dynamics discussed in this paper. The 
basic idea is to modify the ID trapping potential with a 
small but localized attractive potential, which creates a 
tight potential minimum along the weak confining axis 
of the ID trap. 

To create the situation described in this paper one 
needs exceptionally smooth trapping potentials as ob- 
tained for example with atom chips nano-fabricated in 
gold layers on semiconductor substrates [lol l36l IBH l52l| . 
On these atom chips, we can routinely create ID traps 
with an aspect ratio of larger than 1000 53] . The lon- 
gitudinal confinement in the traps can be harmonic or 
box-like on the energy scale of the chemical potential, 
depending on the actual setup. 



1. Creating the dimple with magnetic fields 

On an atom chip, the ID trap is formed in a magnetic 
minimum created by superimposing the magnetic field of 
a current in a straight wire by a homogeneous bias field 
Hi! EH El We assume the trapping wire to be paral- 
lel to the 2-direction (Fig. and the homogeneous bias 
field -Bi-,j as = (£? x , By, -Bj g- e ) to be nearly orthogonal, 



D 



Ioffe 



< | B ± | with B_l = (S x ,S y ,0). The 



i.e. 

component of the field along the wire direction, the Ioffe 
field -Bj ff e , defines the minimum of the trapping field. 

A dimple can be created by a second wire, located 
along the ^-direction, and thus crossing the trapping wire 
(Fig. 56| . A small current jpj in this dimple-wire cre- 
ates a magnetic field which either subtracts from, or adds 
to, -Bi ff e , thus creating a dimple or a barrier respec- 
tively. For a linear current crossing the trapping wire at 
z = 0, with the ID trap at height h above it (Fig. 0(a)), 
the dimple potential is given by 



V D (z) = V 



1 + {z/hf 



(Al) 



The depth of the dimple potential Vq = 
3FMs m i r (A t o/27r)(jr)//i) is determined by the mag- 
nitude of the magnetic field created by the current jrj at 
the position of the trap. Here /is is the Bohr magneton, 
gp the Lande factor of the trapped atomic state \F,mp) 
and /xo = 47r Gmm/A is the vacuum permeability. The 
longitudinal size (i.e. width) of the dimple is set by the 
distance h. 

The geometry of the dimple depends on the ratio be- 
tween h, the height of the trap above the dimple wire, 
and c?trap' ^ ne distance of the atoms to the trapping wire: 



B 



-TJ d trap V 2B Ioffc S D 



(A2) 



The parameters h and dtrap can be controlled indepen- 
dently, by either mounting the dimple wire above the 
trapping wire (Fig. 0a)) or by rotating the trap towards 
the surface by rotating the bias field (Fig. 0b)), thus 
yielding complete control of the dimple trap anisotropy. 
When the ID gas is trapped relatively far from the 
trapping wire, the remaining residual disorder potentials 
lid - l3fij | are reduced to a point where they are not 
visible even for tight transverse confinement of typical 
frequency lu ± ~ 27r(5 kHz) jHH^I- 

In a real experimental setup, the width of the wires and 
the exact form of the current distribution in the wires 
determines the potential. By measuring the magnetic 
dimple potential, as described in Q, we found that it is 
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FIG. 7: (color online) Wire configurations for creating a dim- 
ple in a quasi-lD gas on an atom chip, (a) The dimple is 
created by the crossing wire placed above the trapping wire. 
The magnetic field created by the current flowing in the cross- 
ing dimple wire subtracts from the minimum of the ID trap 
(Ioffe field Bjoffg). (b) By rotating the bias field the ID trap 
is displaced and the dimple is created on the side of the trap- 
ping wire. To keep the trap symmetric around the dimple, 
a wire crossing is again required, (c) In the rotated geome- 
try charging an electrode creates a strongly confining dimple, 
without the need for wire crossing. 



by suitably choosing the position of the trap, as shown 
in Fig. EJc). 

Such a ID electric dimple can be easily created on the 
atom chip with lu±_ ~ 2tt(5 kHz) described above, once 
the trap is rotated towards the chip surface and posi- 
tioned above an orthogonal wire. In the electric case, a 
few Volts of electric potential on the dimple wire create 
an adjustable dimple, whose parameters are similar to 
the ones described in this paper for typical trap heights 
of order h ~ (5 — 25)/xm [57j . 

In a real setup one must again take into account the 
detailed geometry, the exact boundary conditions and the 
dielectric properties of the chip substrate. Measuring the 
potentials we found that they can be described, to good 
accuracy, by a gaussian shape |57j . 

Creating the dimple with electric fields has the ad- 
vantage that a switchable dimple can be created even 
without the need for crossing the wires on the chip, as 
shown in figure E{c). In addition, high resistivity ma- 
terials can be used for charged structures, and this will 
dramatically reduce the thermal noise induced by John- 
son noise currents [HE EE 0, EE EE ESI EE EE E3 > sucn 
that distances of 1 /mi to the charged dimple wire are 
feasible. 



actually very well described by a gaussian shape, as used 
in the calculations presented in this paper |57j . 

On an atom chip, a ID trap with uj±_ ~ 2ir(5 kHz) can 
be typically created by a 1000 mA current and a B± ~ 
28G. If the trap is then rotated towards the chip surface 
and positioned above an orthogonal wire at heights h ~ 
(5 — 25)/im, a current jrj ~ (10 — 100)fiA in the dimple 
wire creates potentials like the ones described in the main 
part of this paper. The larger values of h correspond to 
wider dimples and require larger j-Q . 



2. Creating the dimple by adding electric fields 

A second way of implementing a controllable dimple 
on an atom chip is by using electric fields (Fig. EJc)). 
An atom with electric polarizability a feels an attractive 
potential V cl — — ^aE 2 . Charging up a wire positioned 
orthogonal underneath a trap will create a dimple poten- 
tial |58|. 

In a simple model of a ID trap at height h above the 
line charge Q, the electric dimple potential 



3. Creating a dimple by superposing a dipole 
potential 

A third way of creating a dimple on a atom chip is 
following j2(J and superimposing a dipole potential on the 
atom chip trap E^|. The dimple is created by a tightly 
focused red-detuned laser beam. The size of the dimple 
is given by the focus size of the beam. The form of the 
dimple potential is gaussian, Vq — Vq e - ^ -20 ' l 2a , and 
the size of the dimple is set by the numerical aperture of 
the imaging lens, and will typically be of order of 3-10 
/j,m. 

The dipole dimple can be created at arbitrary distance 
from the surface of the chip. It is the simplest way to 
create a tight dimple in cold atom experiments, and has 
been used in the condensate growth experiments at MIT 
[1E|, and the 133 Cs experiments at Innsbruck jsEEE]. In 
an atom chip environment, one has to take care of the 
stray light from the chip surface and the resulting inter- 
ference structure inside the dimple. Scattered light will 
create speckles and form a additional disorder potential. 



(A3) 



4. Switching the dimple 



has the same form as the magnetic dimple discussed 
above (Eq. QAip ). The ratio of the longitudinal to trans- 
verse confinement in the dimple depends again on the 
distance c^rap between the trapping wire and the atom 
trap, and on the height h above the line charge. This 
ratio can thus be adjusted over a broad range of values, 



All three examples of dimple traps discussed above 
can be switched very fast. In the electric and magnetic 
cases, only very small currents or charges are needed, and 
rapid switching times (r <C 1/xs), much faster than the 
timescale of the transverse trapping frequency lu±, can 
easily be obtained. The same can be said for the dipole 
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dimple, for which the switching depends only on the light switching time, which can be again r <C l/us. 
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